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ABSTRACT 

We study the alignment of grains subject to both radiative torques and pinwheel torques while 
accounting for thermal flipping of grains. By pinwheel torques we refer to all systematic torques that 
are fixed in grain body axes, including the radiative torques arising from scattering and absorption of 
isotropic radiation. We discuss new types of pinwheel torques, which are systematic torques arising 
from infrared emission and torques arising from the interaction of grains with ions and electrons in 
hot plasma. We show that both types of torques are long-lived, i.e. may exist longer than gaseous 
damping time. We compare these torques with the torques introduced by E. Purcell, namely, torques 
due to II2 formation, the variation of accommodation coefficient for gaseous collisions and photoelec- 
tric emission. Furthermore, we revise the Lazarian & Draine model for grain thermal flipping. We 
calculate mean flipping timescale induced by Barnett and nuclear relaxation for both paramagnetic 
and superparamagnetic grains, in the presence of stochastic torques associated with pinwheel torques, 
e.g. the stochastic torques arising from II2 formation, and gas bombardment. We show that the com- 
bined effect of internal relaxation and stochastic torques can result in fast flipping for sufficiently small 
grains and, because of this, they get thermally trapped, i.e. rotate thermally in spite of the presence 
of pinwheel torques. For sufficiently large grains, we show that the pinwheel torques can increase the 
degree of grain alignment achievable with the radiative torques by increasing the magnitude of the 
angular momentum of low attractor points and/or by driving grains to new high attractor points. 
Subject headings: Polarization -dust extinction -ISM: magnetic fields 
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1. INTRODUCTION 

Polarization of radiation arising from emission or ab- 
sorption by aligned grains is widely used to study mag- 
netic fields in the diffuse interstellar medium (see Good- 
man et al. 1995), molecular clouds (see Hildebrand et 
al. 2000, 2002) and in prestellar cores (Crutcher et al. 
2004) . The grain alignment involves the alignment of the 
axis of major inertia sl± with angular momentum J and 
the alignment of J with respect to the ambient magnetic 
field. 

An understanding of grain dynamics is extremely im- 
portant for both understanding of dust properties and 
dust alignment. The latter is essential for determining 
situations when polarization of starlight passing through 
dusty interstellar g well as polarization of dust emis- 
sion can be reliably interpreted in terms of magnetic 
field direction. In terms of dust properties, fast rota- 
tion should disrupt loose aggregates, placing limits for 
the fractal dimensions of dust particles. 

As dust scatters or absorbs photons, it experiences 
radiative torques. These torques can be stochastic or 
systematic. Stochastic radiative torques arise from, for 
instance, a spheroidal grain randomly emitting or ab- 
sorbing photons. The latter process, for instance, was 
invoked by Harwit (1970) in his model of grain align- 
ment based on grains being preferentially spun up in the 
direction perpendicular to the photon beam. However, 
Purcell & Spitzer (1971) showed that randomization aris- 
ing from the same grain emitting thermal photons makes 
the achievable degree of grain alignment negligible. In 
fact, they showed that random radiative torques arising 
from grain thermal emission is an important process of 
grain randomization. More recently, grain thermal emis- 
sion was analyzed as a source of the excitation of grain 



rotation as well as its damping in relation to the rotation 
of tiny spinning grains that are likely to be responsible 
for the so-called anomalous foreground emission (Draine 
& Lazarian 1998). 

Systematic radiative torques were first introduced by 
Dolginov (1972) in terms of chiral, e.g. hypothetical 
quartz grains. Starlight passing through such grains 
would spin them up. Later, Dolginov & Mytrophanov 
(1976) considered an irregular grain model, consisting of 
two twisted spheroids, made of more accepted materials, 
e.g. silicate grains, and claimed that these grains will 
be both spun up and aligned by radiative torques arising 
from the anisotropic component of radiation field. ^ Fur- 
ther on, following the convention we adopted in Lazarian 
& Hoang (2007a, hereafter LH07a), we shall call these 
torques RATs. The work by Dolginov & Mytraphanov 
(1976) was, unfortunately, mostly ignored for 20 years. 
A possible explanation of this may be due to the fact, 
that for a long time, the magnitude of RATs for realistic 
irregular grains remained unclear. 

A renewed interest to RATs was induced by the pos- 
sibility of calculating them for arbitrary grain shapes. 
This occurred after Bruce Draine modified correspond- 
ingly his publicly-available Discrete Dipole Approxima- 
tion code (hereafter DDSCAT; Draine & Flatau 2004). 
Moreover, Draine & Weingartner (1996, 1997, henceforth 
DW96, DW97, respectively) conjectured that RATs may 

^ Using so-called Rayleigh-Gans approximation, Dolginov & 
Silantiev (1976) provided calculations of the radiative torques for 
a model of an irregular grain consisting of two ellipsoids twisted 
with respect to each other, but our calculations of radiative torques 
using the Discrete Dipole Approximation code (Draine & Flatau 
2004) for the same set parameters and the same model as in the 
aforementioned work are in conflict with their analytical findings 
(see more in Hoang & Lazarian 2008b). 
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provide the primary alignment mechanism for interstel- 
lar grains. Support for this claim came through later 
research in the field, by better understanding the dy- 
namics of grains subject to RATs from anisotropic radia- 
tion fields, calculating grain alignment with an analytical 
model (LH07a), as well as including important physical 
processes like grain wobbling (Lazarian 1994; Lazarian 
& Roberge 1997; Weingartner & Draine 2003; Hoang & 
Lazarian 2008a), and accounting within the model for 
gaseous bombardment and uncompensated torques aris- 
ing from H2 formation, as was done in Hoang & Lazarian 
(2008a, henceforth HL08a). As it was described in a re- 
view by Lazarian (2007) , RAT alignment mechanism has 
become not only the leading candidate to explain inter- 
stellar grain alignment, but also to explain grain align- 
ment in many other astrophysical environments, includ- 
ing circumstellar regions (Aitken et al. 2002), accretion 
disks (Cho & Lazarian 2007), comet atmospheres (see 
Roscnbush et al 2007), and molecular clouds (Whittet et 
al. 2008). 

In LH07a, we subjected to scrutiny the properties of 
RATs. Using a simple analytical model (AMO) of a heli- 
cal grain we studied the properties of RATs and the RAT 
alignment. The results obtained by the AMO were shown 
to be in good correspondence with numerical calculations 
obtained by DDSCAT for irregular grains. Invoking the 
generic properties of the RAT components, we explained 
the RAT alignment of grains in both the absence and 
presence of magnetic fields. Intentionally, for the sake 
of simplicity, in LH07a we studied a simplified dynami- 
cal model to demonstrate the effect of RATs. Within the 
latter model we showed that RATs can align grains at at- 
tractor points with low angular momentmn (i.e., J of the 
order of thermal angular momentum Jih, hereafter, low- 
J attractor points), and/or attractor points with high 
angular momentum (i.e., J » Jth, hereafter high- J at- 
tractor points). The high- J attractor points mostly cor- 
respond to a perfect alignment of J with respect to B, 
while the low-J attractor points occur at perfect align- 
ment angle or at some angle in the vicinity of the perfect 
alignment. 

One of the effects that was not discussed in LH07a was 
the effect of thermal fluctuations. Lazarian (1994) no- 
ticed that in spite of the fast rates of internal relaxation, 
rotating grains wobble. This conclusion was a conse- 
quence of the Fluctuation-Dissipation Theorem (see Lan- 
dau & Lifshitz 1976), which states that any dissipation 
process, e.g. internal relaxation, should be accompanied 
by the proportionally fluctuation process. It is possible 
to show that grains rotating thermally are expected to 
wobble with larger amplitude compared to their coun- 
terparts rotating at suprathermal (much greater than 
thermal) rates (see Lazarian 1994, Lazarian & Roberge 
1997). Therefore the zero value of angular momentum 
at low-J attractor points obtained in LH07a stemmed 
from the assumption of the perfect coupling of grain axis 
of the maximal moment of inertia ai with J, which is 
violated for sufficiently low J. In fact, in HL08a, we 
treated fully the dynamics of the grain alignment by tak- 
ing into account thermal fiuctuations (see also Weingart- 
ner & Draine 2003), and showed that thermal fluctu- 
ations within irregular grains can increase the angular 
momentum of the low-J attractor points from J = to 
J Jd, i.e. the angular momentum of grains corre- 



sponding to dust grain temperature. In addition, HL08a 
showed that gas bombardment and other randomizing 
torques can contribute significantly to the RAT align- 
ment by moving grains from low-J to high-J attractor 
points, whenever the latter are present. However, we 
feel that a comprehensive study on the degree of align- 
ment as a function of grain size and radiation intensity 
is very essential for polarization modeling. Note that 
earlier papers which attempted to introduce grain align- 
ment into polarization modeling (Cho & Lazarian 2005; 
Pelkonen et al. 2007; Bethell et al. 2007) assumed per- 
fect alignment for grains larger than some critical size, 
and no alignment for smaller grains. The latter was in- 
ferred using the criterion of whether RATs for radiation 
field under interest, which are sharp function of grain 
size, can induce grains to rotate several times faster than 
the thermal rotation rate. This is a rather crude criterion 
because it does not take into account the type of RAT 
alignment, i.e with or without higli-J attractor points. 
For instance, if RATs align all grains at low-J attractor 
points, then assumption above leads to overestimates for 
the degree of alignment. 

In addition to RATs, other uncompensated torques act 
on grains (Purcell 1979; LD99a; Roberge & Ford 2000, 
henceforth RFOO). Purcell (1979) proposed three pro- 
cesses that can produce uncompensated torques: hydro- 
gen formation, photoelectric effect and the variation of 
accommodation coefficient. These pinwhecl torques to- 
gether with the paramagnetic dissipation were thought 
to be the major mechanism leading to the alignment of 
J with B. However, the paramagnetic dissipation was 
found to be unable to explain the alignment of param- 
agnetic grains in the weak interstellar magnetic field. In 
any case, in the presence of RATs, the alignment arising 
from paramagnetic dissipation for ordinary paramagnetic 
grains is negligible compared to that arising from RATs^. 

The efficiencies of the Purcell torques decrease if grains 
wobble thermally as we discussed above. The effect of 
thermal wobbling was quantified in Lazarian & Roberge 
(1997) and the results of this study were used by Lazar- 
ian & Draine (1999ab, henceforth LD99ab) to predict 
new effects of grain dynamics, namely, thermal fiipping 
and thermal trapping. The thermal flipping is the effect 
of occasional increase of the wobbling angle beyond 90 
degrees. When this occurs, the torques that are flxed 
in grain body coordinates change their direction. In 
the case thermal flipping occurs fast enough, the Purcell 
torques get averaged out and the grain rotates thermally 

^ Unfortunately, the misconception that RATs act as proxies of 
Purccll's torques and the aUgnment arises from paramagnetic dis- 
sipation is well entrenched in the grain alignment Uterature. This 
claim is based on the DW96 study and, in spite of the fact, that 
the efficient RAT alignment was demonstrated already in DW97 
and later works, the claim persisted. This may be due to the 
fact that, DW97 did not provide, unlike DW96, the predictions 
for the degree of alignment. The problem of predicting the de- 
gree of alignment happened to be a tough one. The present paper 
is an attempt in this direction. Interestingly enough, Dolginov 
& Mytrophanov (1976) did consider RATs as a mechanism capa- 
ble of aligning grains irrespectively of the presence or absence of 
paramagnetic relaxation. For superparamagnetic grains with suf- 
ficient number of iron atom per cluster, Lazarian & Hoang (2008) 
showed that the alignment always happen with high-J attractor 
points and is perfect. However, even in this case, RATs are more 
important than the enhanced paramagnetic relaxation. The latter 
effect causes the stabilization of the high-J attractor point created 
by RATs. 
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in spite of the presence of the uncompensated pinwheel 
torques, i.e it is thermally trapped. 

However, this picture was challenged recently by Wein- 
gartner (2008). His study has several new mathematical 
points. The most important one is that he used a di- 
mensionless variable q — 2IiE/J^ with Ii being the in- 
ertia moment along the axis of major inertia ai and E 
being the total rotational energy to describe the inter- 
nal relaxation, instead of the angle 9 between ai and J. 
He suggested an integration constant for diffusion coeffi- 
cient of internal relaxation, which differs from that in the 
Lazarian & Roberge (1997). We agree with the choice of 
the integration constant. For this choice in Weingartner 
(2008) the diffusion coefficient vanishes when ai perpen- 
dicular to angular momentum J. As a result, he found 
that grains do not experience thermal flipping as a result 
of internal relaxation. We agree with this choice of the 
integration constant as it corresponds to the nature of the 
Fluctuation-Dissipation Theorem employed in Lazarian 
(1994) to describe the wobbling of the grains. Indeed, ac- 
cording to the theorem no dissipation should correspond 
to no fluctuations. For an oblate grain the internal dis- 
sipation goes to zero as 9 Tr/2. Thus the Fluctuation- 
Dissipation Theorem suggests that the at this point the 
fluctuations should also go to zero. One can check that 
the choice of the constant in Weingartner (2008), which 
was obtained from more formal considerations, corre- 
sponds to this physical requirement. However, as it clear 
from the rest of the paper, we disagree that this modifi- 
cation of the diffusion coefficients will preclude physical 
processes of thermal flipping and trapping introduced in 
Lazarian & Draine (1999) from happening. 

In what follows, we consider a more realistic picture of 
suprathermal rotation, namely, we take into account the 
fact that any pinwheel torque is accompanied by stochas- 
tic torques of the same nature. For instance, apart from 
systematic torques, the process of H2 formation induces 
stochastic torques. In this model it is clear that the dif- 
fusion is present even when the angle between ai and 
J is 90 degrees, which means that the thermal flipping 
is possible. In this paper we revisit the LD99a study 
and confirm that the thermal flipping and thermal trap- 
ping predicted in LD99a are important effects of grain 
dynamics. 

The structure of the present paper is as foUowings. 
In §2 we present pinwheel torques, including long-lived 
and short-lived torques. We calculate the maximal value 
of angular momentum that these pinwheel torques can 
achieve for the ISM. In §3 we revise the problem of ther- 
mal hipping and calculate the thermal mean flipping time 
induced by internal relaxation in the presence of exter- 
nal stochastic torques. Critical size of flipping grain and 
trapping size for the ISM are estimated in this section. 
We study the effects of H2 pinwheel torques and thermal 
flipping by stochastic torques accompanied with H2 for- 
mation on the RAT alignment for the diffuse ISM in §4. 
We also calculate the value of J at low attractor points 
as a function of the magnitude of H2 pinwheel torques 
for different grain sizes in this section. We summarize 
our findings in §5. 

2. PINWHEEL TORQUES: REVISITING PURCELL 1979 

Pinwheel torques were first discussed by Purcell (1979), 
who identified H2 torques as the dominant pinwheel 



torques acting in interstellar gas. With the existing 
uncertainties in calculation of the torques this may or 
may not be true. Moreover, for other environments 
(e.g. molecular clouds), we know for certain that it is 
not true. This induces our extended discussion of the 
pinwheel torques. Purcell also estimated that the H2 
torques should be short-lived and this induced further 
researchers (see Spitzer & McGlynn 1979; Lazarian 1995; 
Lazarian & Draine 1997) to concentrate on the studies 
of short-lived pinwheel torques. 

Usually, wc represent the magnitude of pinwheel 
torques in terms of the thermal angular momentum, and 
the gas damping time. For simplicity, we consider a 
spheroid grain, with the moments of inertia /|| — Ii, I± = 
I2 = h-^ They are defined by 



Jth — \J I\\kBTgas — Y ksTgas, 

= 5.89 X 10-20a5_/25i/2^i/2f V2 g cm^rad s'\ (1) 

with p density of material within the grain and p — 
p/3 gcm~^, a grain size, a_5 = a/10~^ cm, s = s/0.5 

with s = b/a and Tgas gas temperature and Tgas — 
Tgas/WQK. And the damping time due to gas colli- 
sion (see Roberge, DeGraff & Flatherty 1993, hereafter 
RDF93) is 

t^^^^ h , 

® Ay/n nnmHa'^Vth^n' 

= 2.3xl0i2^_5.ir-,y2(^^^^) s, (2) 

where vth is the thermal velocity of hydrogen atom with 
density hh, ttih is the hydrogen mass, Fji is the eccen- 
tricity parameter which is unity for sphere (see RDF93), 
and we adopted standard parameters in the ISM. 

We also define dimensionless angular momentum and 
torque 

J' -7^,, (3) 

Jth 

r' = r%^, (4) 

Jth 

where Jth and tgas are given by equations (H]) and 

2.1. Long-lived pinwheel torques 

2.1.1. Radiative torques due to absorption and scattering of 
isotropic starlight 

While Dolginov & Mytrophanov (1976) considered 
only torques arising from anisotropic radiative flows, 
DW96 noticed that radiative torques are not zero even 
if the radiation field is isotropic. Such "isotropic ra- 
diative" torques have all the properties of the pinwheel 
torques introduced by Purcell (1979), but with the ex- 
ception that these torques were definitely long-lived. In- 
deed, to change them, one required to change the entire 
grain shape, rather than to change properties of grain 

^ For oblate spheroid with major and minor axes a and b, 
7|| = 2Ma^/5 = {STTp/15)ba*, and = (4.71 p/ 15) a^b{a'^ + b'^), 
but sometimes we define the equivalent sphere with the same vol- 
ume of the grain, and use Ii = (87r/15)pa^. 
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surface. Tlie former, in most cases required time-scale 
longer than the gaseous damping time tgas- 

The isotropic radiative torques for studied grain shapes 
are usually a factor of 10^ weaker than the torques aris- 
ing from anisotropic radiation. Therefore, in some sub- 
sequent publications (e.g. in Cho & Lazarian 2005, 2007; 
Pelkonen et al. 2007), they were disregarded. They may 
be occasionally important for grain spin-up. In this pa- 
per, we shall consider isotropic radiative torques as an 
example of long-lived pinwheel torques. 

The isotropic radiative torque is defined by 



Q^so = 4^ /o Jo Qa, (O, /3)sin Qd^dQ, 



(5) 



where O is the angle between the axis of major inertia ai 
and the radiation direction k, (3 is the rotation angle of 
the grain along ai (see Fig. [5]), and (5ai(6,/3) = Qr-ai 
with Qr being the radiative torque efhciency vector (see 
Appendix C) is the component of RATs along ai. Here 
(5ai(0, P) does not depend on the precession angle <I> of 
the grain axis ai about k. 

2.1.2. Radiative torques arising from grain infrared emission 

All systematic* radiative torques so far were assumed 
to arise from absorption and scattering of radiation by 
grain (Dolginov & Mytrophanov 1976, DW96). However, 
it is possible to show that pinwheel torques should also 
arise as a result of grain infrared emission. 

To calculate the radiative torques arising from grain in- 
frared emission, let us consider a grain at thermal equilib- 
rium. In such conditions the grain should not be subject 
to systematic spin-up. This means that the torques aris- 
ing from the infrared emission should compensate 
the torques Fe^t arising from absorption and scattering 
of radiation, i.e. 

Fem -I- Text = 0. (6) 

Therefore, when Fe^.^ is known, equation ^ allows us to 
find the torques arising from grain emission. The pin- 
wheel torque due to infrared emission is then 



Qiso i'-^d) 



(7) 



where a is the grain size, X{Td) and Urad are the mean 
wavelength and energy density of infrared emission from 
the dust with temperature Td, and Qiso is the isotropic 
torque (see Appendix C) averaged over the spectrum of 
grain infrared emission. Here the infrared emission spec- 
trum of dust grain with temperature Td is assumed to 
be the blackbody with the energy density per unit of 
wavelength uy. 

u. = (8) 
c 

where B\{Td) is the brightness of the black body with 
temperature Td and c is the light speed. 

We can calculate the maximum value of angular mo- 
mentum induced by Fem as 



1.2 fimj \uisRF 



Urad \ ^1/2 1/2 



/3000 cm-^K 
V nnTgas 



so{Td) 



(9) 



Harwit (1970), Purcell &: Spitzer (1971) considered stochas- 
tic torques arising from emission of photons (see a more detailed 
discussion in Draine & Lazarian 1999). 
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Fig. 1. — J^?i^x/Jth excited by pinwheel torques duo to infrared 
emission for a grain size a = 0.2/im in reflection nebulae (RN), 
photodissociation region (PDR) and HII regions with dust tem- 



perature = 40, 50 140 K, gas temperature Tg, 



100, 300 and 



8 X lO^K, gas density rtg 



10^,10^ and 5 X 10^ cm ^, respec- 



tively. Astronomical silicate grains of shape 1, 4 and 5 in LH07a 
are adopted. 

where ujsjip is the mean energy density of the interstel- 
lar radiation field (hereafter, the ISRF; see Table 1). 

To obtain OZjTd), we first use DDSCAT (Draine & 
Flatau 2004) to calculate RAT components Qei, Qe2 and 
Qe3, and then Qiso for different angles O ranging from 
to TT and /? in the range to 27r at $ = and for the 
entire spectrum of radiation field with A from 0.1/im to 
20/im, for astronomical silicate grains of irregular shapes 
and of different sizes. Then, we average over the angles 
f] and Q following equation ([5]) to get Qiso- Finally, we 
average Qiso over the entire spectrum of grain emission to 
obtain QisoiTd). Combining with equation ([9]) we obtain 

JmTx- Figure [1] shows J^^x/^th for grain shapes 1, 4 
and 5 in LH07a with the size a — 0.2/im, embedded in 
reflection nebulae (RN), photodissociation (PDR) and 
HII regions with dust temperature Td — 40, 50 and 140 
K, gas temperature Tgas = 100, 300 and 8 x lO^K, gas 
density rigas = 10'^, 10^ and 5 x 10"^ cm""^, respectively. 

Figure[T]shows that the infrared emission torques are of 
marginal importance for RN and PDR when Td is low. 
However, in H II regions, where the dust temperature 
can achieve ~ 140K as a result of the heating from both 
the radiation of O and B stars, 0.2/im grain can rotate 
suprathermally, i.e., with J > Jth (see Fig. [T]for shape 
4 and 5). Since radiative torques increases fast with the 
grain size (CL07, HL07a), larger grains are expected to 
rotate faster by infrared pinwheel torques. Therefore, we 
believe that this important source of pinwheel torques 
should be included in modeling of grain alignment. 

2.1.3. Regular Torques from Grain- Electron interactions in 

plasma 

Heating of grains in plasma is well described in the lit- 
erature (see Dwek et al. 2007). Electrons hit grains more 
frequently than ions. This effect is partially compensated 
by grains getting negative charge. In situations when 
grain heating is mostly by electrons, the heating rate de- 
pends on whether grains are large enough for electrons 
to stop within grain material or electrons pass through 
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grains. Below, we discuss, however, a different effect. 

As electrons pass through grain material they exert 
torques. These torques can have random and systematic 
components. Similarly, as for radiative torques, for ho- 
mogeneous grains with symmetric shapes, e.g. ellipsoids, 
prisms, systematic torques vanish. However, if grains are 
inhomogeneous or of irregular shapes, we expect to ob- 
serve Purcell-type torques arising from grain-plasma in- 
teractions. For instance, electrons penetrating from dif- 
ferent sites of a grain having a shape of prisms, but made 
of inhomogeneous material may have different stopping 
distance, deflection rates or probabilities of ejecting sec- 
ondary electrons. Similar as in the Purcell (1979) study, 
these variations at different sites are difficult to quantify 
with a precise model. Therefore, as in Purcell (1979) one 
may introduce variations of grain properties using plau- 
sibility arguments. One quantitative estimate is beyond 
the scope of this paper. 

2.1.4. Regular Torques due to gas-grain interaction 

Below in §2.2.2 we discuss systematic torques arising 
from the variations of the accommodation coefficient. 
However, gaseous bombardment can produce torques re- 
lated to grain shape. In Lazarian & Hoang (2007b) we 
discussed torques arising from grain streaming^, which 
we assumed helical and approximated with AMO. The 
torques, arising from such a gas-grain interaction are 
analogous to RATs. It is only natural to assume the 
existence of the analog of the isotropic torques for grain- 
gas interactions. In terms of AMO this would correspond 
to the opposite surfaces of the mirror interacting differ- 
ently with the impinging atoms, i.e. very similar to the 
modification of AMO for representing of isotropic radia- 
tive torques that we discussed in LH07a. Such torques 
should be long-lived, as the changes of the grain shape 
would be involved. 

2.2. Short-lived pinwheel torques 

Below, following Purcell (1979) we identify the torques 
due to H2 formation and the variation of the accommoda- 
tion coefficient over grain surface as short-lived pinwheel 
torques. It was argued in Lazarian (1995) that desorp- 
tion of atoms can make these torques long-lived. Never- 
theless, conservatively, we follow Purcell's approach. 

2.2.1. Regular torques arising from H2 formation 

H2 formation followed by the ejection of H2 molecules 
results in the systematic torque along the preferred ro- 
tation axis (i.e., axis of major inertia ai). In addition, 
the H2 formation on different active sites is random, and 
leads to the second order effect of torques. 

The detail calculations for the systematic torque com- 
ponent for a brick and a spheroid were presented in Pur- 
cell (1979) and Lazarian & Roberge (1997), respectively. 
For a brick grain of sides 2a x 2a x 26 with s = b/a < 1, 
the H2 pinwheel torque in dimcnsionlcss units is given 

by 

((r'^J^)V2«3.5x 10^(1 -y)7^J^/^ai4^pV2f-;/^ 

Unlike the traditional mechanical alignment processes (Gold 
1952b, Purcell 1969, Purcell & Spitzer 1971, Dolginov & Mytro- 
phanov 1976, Lazarian 1994, 1995b, 1997, Lazarian & Efroimsky 
1996, Efroimsky 2003, Roberge et al. 1995), this process is efficient 
even if the gas-grain drift is subsonic. 



x^.»(^^j ' (10) 

where y — 2n{H2)/nH being the fraction of hydrogen 
in the molecular state, jh = 7i//0.1 with being the 
fraction of H atom converted to H2 molecule, s = s/0.5, 
a is the density of active site on the grain surface and 
Ekin = £'fci„/0.2eV with Ekin being the kinetic energy 
of H2 molecule escaping from the grain surface. We can 
now write the maximal value of grain angular momentum 
achieved by H2 pinwheel torque as 

- {{t'hS)'"- (11) 

Equation (fTU|) shows that the H2 pinwheel torque de- 
pends on the active site density a. This parameter is 
not well constrained at the moment. Lazarian (1995) 
argued that poisoning of active sites is more important 
than resurfacing and estimated that small grains with 
a < 10~^ cm should have active sites poisoned. 

2.2.2. Regular torques due to the variation of 
accommodation coefficient 

Purcell (1979) noticed that if gas and grain temper- 
atures are different, the variation of the properties of 
the dust grain should create systematic torques. For the 
variation of the accommodation coefficient of 10 percent 
with the scale of 1/3 of the grain size Purcell (1979) ob- 
tained that for typical conditions of the diffuse medium 
the torques will be of the order 10~^ of those by H2 
torques. 

In molecular clouds where torques due to H2 formation 
are inefficient due to the lack of atomic hydrogen, grains 
can rotate suprathermally by the torques discussed in 
this section, or by isotropic radiative torques. This agrees 
with observations that indicate that in molecular clouds 
the temperatures of dust and gas differ (see Lazarian, 
Goodman & Myers 1997). 

2.3. Comparison to Anisotropic Radiative Torques 
2.3.1. Anisotropic Radiative Torques 

An analytical model (AMO) of RATs arising from the 
perfect refiection of the radiation beam onto the sur- 
face of a helical grain was proposed in LH07a. This 
grain model consists of a perfect reflecting mirror, con- 
nected to an ellipsoid (see the upper panel in Figure [2]). 
LH07a proved that the AMO represents very well the ba- 
sic properties of RATs for irregular grains obtained from 
DDSCAT. 

Throughout the paper, we adopt the functional forms 
of RATs from the AMO, which are given in Appendix 
C. To study grain alignment, the amplitude of torque 
components for the AMO, described by g""*^^ can be 
rescaled to the values for irregular grain obtained us- 
ing DDSCAT. To account for the magnitude of RATs as 
functions of grain size and wavelength of radiation, we 
make use of the approximate scaling obtained in LH07a, 
which is given by 

Qr«0.4Qy, (12) 

where 7; = for X ^ a and rj = —3 for A 3> a (see 
Appendix C for more detail). 
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Fig. 2. — Upper panel: Geometry for AMO consisting of a per- 
fectly reflecting spheroid and a weightless square mirror of side I2 
tilted by an angle i, connected to the spheroid by a rod with length 
'2 S> ^2- Lower panel: The orientation of a grain, described by 
three principal axes ai , 02 , 03 , in the laboratory coordinate system 
(scattering reference system) 61,62,63 is defined by three angles 
0,/3, The direction of incident photon beam k is along e\. 



The magnitude of RATs in dimensionless units is 

YRAT_j^i _ l^Urada^tgasyr- „x 

J^^T — i R4T — T—: ^r^I \^'^) 



^2 X 10* 



A 



.2 



UISRF 



'fJ 



flHTn. 



(14) 



where 7 is the anisotropy of radiation field, A and Qp 
are the wavelength and RAT efficiency averaged over the 
spectrum of radiation field. 
Using equation for Qp by replacing A = A, we get 



J, 



'RAT 



.8 X 10^ f^) ( ^—T' ( ^) pV2„i 
Wl)\1.2 ^lm) \uisrf) ' 

(15) 



/2 

-5 




Fig. 3. — Alignment coordinate system: i/i is the angle between 
the magnetic field B and the radiation direction k, § is the angle 
between the angular momentum vector J and B, and cf> is the Lar- 
mor precession angle of the magnetic moment parallel to J around 
B. 

2.4. H2 formation torques versus anisotropic RATs 

To study the relative role of pinwheel torques® to 
RATs, let us define r^f2/r_RAT to be the ratio of magni- 
tude of H2 torques to RATs. Using equations (fTO)) and 
(fT^l). we obtain 



j'H2 



r' 

max H2 



PRAT p' 







RAT 



ri+1 



^ ^0.1\ ^a-5y ^ l■2^JLm \ ■ uisRF 



X nHf^l^s'^'{l-y)jH 



IQi 



1/2 



E 



rad 

1/2 

kin ' 

(16) 



The H2 torque-RATs ratio depends on both the grain 
size a and the active site density a, which is a rather 
uncertain parameter. Therefore, we consider a range of 
a with the lower limit amin ~ l/47ra^, corresponding to 
one active site on the entire grain surface, and the upper 
limit a = 10^^ cm~^, which is a reasonable value for high 
density of active sites. For a_5 = 1 grain, equation (|16p 

gives J'maxl J'm^ 10^ and for a = lO^^cm"^ using 

standard parameters of the ISM. 

Let us calculate the ratio of maximal angular mo- 
mentum to its thermal value for the H2 formation (eq. 
[TU|) . and for the anisotropic RATs Qr calculated using 
DDSCAT for a more realistic grain (shape 1). Qp is ob- 
tained by averaging Qp over the spectrum of ISRF as 
in Mathis, Mezger & Panagia (1983). Figure |4] shows 
JmaxI Jth for these two torques. Here we assume that 
radiation is parallel to the magnetic field and disregard 
the decrease of RATs due to thermal wobbling of the 
grain. The shaded area shows Jmax/Jth for H2 torques 
where the upper and lower boundaries correspond to the 
lower and upper limits of density of active sites a. 

It can be seen that Jmax/Jth for RATs increases with 
grain size faster than that due to H2 pinwheel torques, 

® In the rest of the paper, we consider H2 torques as a represen- 
tative of pinwheel torques. 
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Fig. 4. — The maximum value of angular momentum relative to 
the thermal value that H2 pinwheel torques and RATs can achieve 
as a function of grain size for the ISM. The shaded area repre- 
sents the pinwheel torques due to the H2 formation (eq. IIUI I. and 
the solid line shows represents RATs for an irregular grain (shape 
1) and the ISRF. For H2 pinwheel torques, the upper boundary 
corresponds to one site per a grain of a = 10~®cm, and the lower 
boundary corresponds to an reasonable value for high density of ac- 
tive sites. For RATs, radiation direction is assumed to be parallel 
to the magnetic field. 

and becomes dominant for grains larger than 5 x 10~^cm. 
In contrast, for grains smaller than 7 x lO^^cm, H2 pin- 
wheel torques are dominant over RATs. ^ 

In addition, such small grains can still rotate 
suprathermally with Jmax S> Jth by II2 torques. In 
principle, the suprathermal rotation of small grains can 
produce a detectable degree of alignment as a result of 
paramagnetic dissipation. This contradicts to observa- 
tions that grains smaller than 5 x 10~^cm can not be 
aligned (Kim & Martin 1995). This contradiction can 
be resolved if grains undergo thermal flipping and get 
thermally trapped, as we discuss below. 

We note that the value of RATs for a given grain size, 
depends on the angle "0 between the magnetic field B and 
the radiation direction k. Our study in HL08b shows that 
the value of the spin-up component of RATs at the angle 
ijj approaching 7r/2 may be 10~^ of the value suggested 
by equation (jl6p making RATs weaker than the pinwheel 
torques for some angles tj}. 

3. THERMAL FLIPPING AND TRAPPING: REVISITING 
LAZARIAN & DRAINE 1999 

3.1. Barnett, Nuclear and Superparamagnetic Internal 
Relaxation 

Barnett (1915) pointed out that a rotating paramag- 
netic body gets magnetized with the magnetic moment 
parallel to the angular velocity^. This effect can be easily 
understood using a classical model. 

Consider a paramagnetic grain which is rotating with 

For a_5 = 1 grain, Figure |4] shows that Jmax/ Jmax^ ~ 60 
for a = lO^^cm"'^. This value is much smaller than that given 
by equation I II6I I. Such difference stems from t he fa ct that Qp for 
the ISRF is larger than that given by equation II12I I used to derive 
equation III6I I. 

* This is an inverse of the Einstein-de Haas effect, that was used 
to measure the spin of the electron. To grain alignment theory the 
Barnett effect was introduced by Dolginov & Mytrophanov (1976), 
who noticed that the effect should induce the magnetic moment of 
grains. 



the angular velocity fl. As the grain rotates with angu- 
lar velocity $7, the torque acting on the electron spin is 
J^J = r2 X J. The equivalent torque can be induced by 
a magnetic field Hb acting on the magnetic moment /i 
associated with the spin, i.e. 1/cfj, x Hs, which provides 
the Barnett-equivalent field Heqv = which can also 
be presented as 



19 



(17) 



where 7^ = ^/{Jc) = eJ/(meJc) = e/{raec), with e 
electron charge and me electron mass is the magneto- 
mechanical ratio of an electron. 

Purcell (1979) showed that as the grain wobbles, the 
changes of magnetization in the grain axes cause the in- 
ternal relaxation, which he termed " Barnett relaxation" . 
LD99b revisited the problem by taking into account both 
spin-lattice and spin-spin relaxation (see Morish 1980). 
The resulting rate of change for the angle 9 between the 
grain axis of maximum moment of inertia ai and J (see 
Figure [5]) is 



VK{h-l)f' 



sin dcosdf (9, uj) (18) 



where V is the volume of the grain, h = I\\/I±, K — 
Xe/wi = x(0)Te;/[l + (wiTe/)2]2 being the imaginary 
part of the magnetic susceptibility with uji — Jcos 9 , 

X(0) = 4.2 X IQ-'^fpf-^ with /p being the fraction of 
paramagnetism material in the grain, and Tei is the re- 
laxation time of electronic spins (Draine 1996; LD99a). / 
is the factor of order of unity, and when / = 1 the result 
coincides with that of Purcell (1979). Ignoring the fac- 
tor / in the above equation, one gets the Purcell (1979) 
expression for the Barnett relaxation rate for a brick as 
in Figure[5]with I» = lQ/?,pa%, h = 2aV(a^ + b'^)- 



tsa 



VKh^{h-l)J'^ 
;2x 10Va-5s[0.5 + 0.125sT ( 



(19) 



where Tf.i ~ r2 ~ 2.9 x 10~^^/j7^s with assumption of 
fp = 0.1 is the spin-spin coupling time, p = p/(3 gcm"'^) 
and d_5 = a/10^^ cm; Jd = yJI^ksTd/ [h — 1) is the 
dust thermal angular momentum ^. 

Although Purcell (1979) considered grains with both 
electron and nuclear spins, his study missed the effect of 
internal relaxation related to nuclear spins. LD99b found 
that for astrophysical grains of realistic composition nu- 
clear spins induce a new type of relaxation, which was 
termed by LD99b "nuclear relaxation". This relaxation 
can be understood in a simple-minded approach in terms 
of much stronger equivalent magnetic field given by equa- 
tion (fT7|) . Indeed, this field is proportional to the mass 
of the species involved. As the paramagnetic relaxation 

^ The relaxation of electron spins results from the spin-latice and 
spin-spin relaxation, with time scales ri S> t2, so here we adopted 
Tei ~ T2 (Draine 1996) 
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Fig. 5. — Grain model used for calculations of internal relaxation. 

is proportional to H^q^x" ^ IgXi^)''' ^ m?{l/m^)T ^ t, 
and as t„ ^ Tez, one can understand the nature of the 
dominance of the nuclear relaxation, which has a char- 
acteristic time: 

2 



in 



cl '■ 



6 X lO^p^^Ls 



1 



Jd 
J 

s, 



5[0.5 + 0.125s 



2i2 



(20) 



where r„ ^ = r„„^ + r„g^ is the relaxation rate induced 
by the nuclei-nuclei and electron-nuclei spin interactions 
(see LD99b). 

Starting with the classical study by Jones & Spitzer 
(1968) superparamagnetic inclusions are frequently con- 
sidered as a constituent part of astrophysical dust grains 
(see Mathis 1986; Martin 1994, Goodman et al. 1995; 
Draine & Lazarian 1999; Roberge & Lazarian 1999). 
However, the effect of superparamagnetic inclusions on 
the internal relaxation was considered only recently, 
namely in Lazarian & Hoang (2008, LH08). There it 
was shown that superparamagnetic inclusions enhance 
the efhciency of the Barnett relaxation and also extend 
the range of frequencies for which the nuclear relaxation 
is efficient. In particular, the superparamagnetically en- 
hanced rate of the Barnett relaxation is 



is, 



lO^p^aL^—exp 



1 



ar,sup ' 



-O.OllA^d ( Jd 



J 



2l2 



xs[0.5 + 0.125s^] 



\ 2 



s(21) 



where r"^ = lifey.Y>{-NciTact/Td) with Na being the 
number of iron atom per cluster and Tact — O.OllK is 
the rate of remagnetization due to superparamagnetic 
inclusions, and the fraction fsup — 0.01 of atoms being 
magnetic was adopted (see LH08). 

In addition, Purcell (1979) considered the internal re- 
laxation arising from grains being inelastic. The im- 



FlG. 6. — Timescales for different internal relaxation processes 
as a function of the grain size for J = SJ^j and J^. The number of 
iron atom per cluster N^i = 2 X 10^ is assumed. 

proved rates of this relaxation were obtained in Lazar- 
ian & Efroimsky (1999). However, for the typical grains 
in the diffuse interstellar medium, this rates are sub- 
dominant. The inelastic relaxation gets important for 
large grains, which we consider in other papers (e.g. in 
HL08b). 

All in all, the total internal relaxation rate t,^j(J) is 
the sum of all relaxation rates 



tl 



Bar, sup ' 7iuct 



(22) 



where ... are placed instead of other relaxation rates, e.g. 
inelastic relaxation rate. Figure [5] shows the timescales 
for various internal relaxation processes as a function of 
the grain size for J = 5Jd and Jd- We assume Nd = 10^ 
for superparamagnetic inclusions. 

3.2. Revised model of thermal flipping 

Internal relaxation induces the thermal wobbling 
(Lazarian 1994; Lazarian & Roberge 1997), which ac- 
cording to LD99a can result in interesting effects like 
thermal flipping and thermal trapping. An alterna- 
tive, but consistent with LD99a treatment of the ther- 
mal flipping was introduced in a preprint by Roberge 
& Ford (2000; hereafter RFOO). Both treatments use 
the diffusion coefficients derived in Lazarian & Roberge 
(1997). However, Weingartner (2008) questioned both 
the derivation of the diffusion coefficients and the phe- 
nomena of thermal flipping and thermal trapping. 

Here we introduce a model that generalizes the one 
in LD99a, namely, it takes into account the fact that 
the processes that induce pinwheel torques also induce 
stochastic torques. Such torques should be taken into ac- 
count in the treatment of suprathermally rotating grains. 
In addition, astrophysical grains are embedded in gas, 
that also induces stochastic bombardment. Formally, 
the diffusion coefficients arising from the internal relax- 
ation and external stochastic torques are different. For 
instance, internal relaxation cannot change the angular 
momentum, but external torques can do so. However, 
within the formalism that we are considering, the exter- 
nal stochastic torques are subdominant for all 9 angles 
apart from 6 tt/2. Therefore, for the sake of simplicity, 
we disregard the small change of the angular momentum 
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TABLE 1 

Physical parameters used for calculations in the paper 



Flipping by flue, ISM, J = J 



Parameters 



Diffuse ISM 



MCs 





30 




100 


Td(K) 


15 


y=2n(H2)/nH 


1 


7 


0.1 


A 


1.2 




8.64 X 10" 


7-ff 


0.1 


igas (s) 


2.3 X 10^2{ 



r)a- 



300 

20 

15 

10-3 

0.1 

1.2 

varied 
0.1 

the same 



n/30 g cm ^, 
'3 where p = 3 g cm" 



and 
3 ;c 



Note. — Here Tgas = T, 
a_5 = agff/10~^ cm. p = p/3 g cm 
the mass density of the grain. 

amplitude that the external stochastic torques cause. In- 
stead, we concentrate on their influence on the variation 
of the angle 0.^° 

Combining dimensionless diffusion coefficients for in- 
ternal relaxation, and impulses due to H2 formation and 
gas bombardment, we get 

t„ 



B' 



H2 



B'cou) 



(23) 



where B'{9) and i?^^ are given by equations (|A6P and 
(jA13|) . and B'^As obtained using equations (|A13p com- 
bined with (jA14[) - (|A16|) . The relative role of internal 
relaxation and H2 formation and gas bombardment are 
described by the ratio tgas /tint- 

Our treatment of grain flipping is similar to the treat- 
ment in RFOO^^, where the approach to solve analyti- 
cally the stochastic differential equations was adopted. 
That approach was described in the textbook of Gar- 
diner (1983, §5.2.7). For the sake of completeness, we 
discuss this approach in Appendix B. Gardiner (1983) 
described the problem of the escape of a particle initially 
confined within the range [a, b] with x = a being a reflect- 
ing boundary and x = h being an absorbing boundary. 
In our flipping problem, the escape corresponds to the 
angle 9 passing 7r/2 value. When 9 gets larger than 7r/2, 
following RFOO, we claim the flipping event. 

In the following, we adopt this approach but study the 
effect of thermal flipping in the presence of impulses due 
to II2 formation and gas bombardment. As we discussed 
earlier, the internal relaxation is fast, so it is reasonable 
to assume that angular momentum is constant in the 
presence of external torques. We calculate the mean flip- 
ping timescale for all internal relaxation processes. For 
the superparamagnetic case, the number of iron per clus- 
ter Nci — 2x 10"^ is adopted. Because the mean flipping 

Ideology of this approach is similar to that employed by 
Spitzer & McGlynn (1979) and Lazarian & Draine (1997) while 
treating grain crossovers. 

The treatment in LD99a is based on the Zeldovich approx- 
imation to the first order transitions. It can also be generalized 
to the presence of the external torques, but the generalization is 
somewhat involved. 
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Fig. 7. — Mean flipping time relative to the gas damping time 
of grain induced by internal relaxation in the presence of impulses 
associated with H2 formation for the ISM (upper) and gas bom- 
bardment for molecular cloud (lower). The angular momentum is 
assumed to be constant and J = Jd- Physical parameters used in 
calculations are given in Table 1. 



time depends on the initial angle 9 (see eq. IB17|) . here 
we show the mean flipping time for 9 = 30°. 

When we disregard the external stochastic torques, we 
flnd, in agreement with Weingartner (2008), that the 
mean flipping time is inflnite, i.e., the grain can not flip 
by internal relaxation only. However, for any realistic 
rates of external, e.g. stochastic torques, the flipping is 
present. 

Figure [7] shows the ratio of mean flipping time to the 
gas damping time, ttf/tgas, as a function of the grain 
size a for the diffuse ISM (upper panel) and a molecular 
cloud (MC; lower panel). For the ISM, the stochastic 
torques due to H2 formation are dominant because the 
number density of atomic hydrogen n{H) is comparable 
to the total number density uh-^^ But in the MC, most 
hydrogen is in molecular state, i.e., nn ~ n{H2), so we 
can disregard impulses due to H2 formation and account 
only for gas bombardment. 

For smaller grains, the internal relaxation is much 



We disregard the effect of gas bombardment in 

ISM because of the dominance of H2 impulses, i.e., 

Ekin/'iksTgas ~ 0.2 X 1.610-12/(4 X 1.38 x IQ-") ~ 6 > B,^„, 
1. 



the 
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stronger than the external impulses, Uf/tgas changes 
rapidly with the grain size a, and the curves for Uf/tgas 
for different relaxation processes follow the steep slopes 
as seen with Unt- For larger grains, the internal relax- 
ation gets weaker, and Uf/tgas exhibits slower increases 
with a. 

3.3. Critical size of flipping 

Let us consider first whether the grain undergoing 
spin-down due to pinwheel torques experiences a regular 
crossover in the fashion described in Spitzer & McGlynn 
(1979) and Lazarian & Draine (1997) or it undergoes a 
flip-over. The regular crossover time according to Lazar- 
ian & Draine (1997) is 

tcros — ~7~r TTT — ' ( ) 



Flipping by fluctuations for the IS^ 



dJn / dt 



- pinw 



where Tpinw is the magnitude of pinwheel torques (see 
§2). Here we consider the pinwheel torques due to H2 for- 
mations as a representative for various kinds of pinwheel 
torques discussed in §2. We also assume that J_l — Jd 
as in Lazarian & Draine (1997). Then equation ([M]) can 
be rewritten 

1/2 

,(25) 



tgas J'ynax V '^th 



Jl 
'J ^ 



{h- l)Tgas, 

where J^ax — J max- following, we consider only 

the pinwheel torques due to H2 formation. It is trivial to 
generalize the treatment for other pinwheel torques. The 
crossover time given by equation (|25p will be compared 
with the mean flipping time tt / to get the critical size of 
flipping. 

Using the pinwheel torque due to H2 formation (eq. 
fTO|) with the use of equation ([25]) . we can calculate tcros 
for various grain sizes a and active site density a. For a 
given a, the critical size of flipping, Ocri is then derived 
by setting tcros = ttf, where ttf are calculated in § 3.2. 

The upper panel in Figure [8] shows the mean flipping 
time ttf in the presence of stochastic H2 torques, and 
crossover time tcros for three values of active site den- 
sity a = 10^^, lO^'^ and 10^^ cm~^. The solid line shows 
ttf when all internal relaxation processes (Barnett, nu- 
clear and superparamagnetic inclusions) are accounted 
for, and the dot line represents that when only Barnett 
and nuclear relaxations are taken into account. The filled 
circles denote the intersects of ttf and tcros, i-c where 
ttf = tcros- The grain size corresponding to this intersect 
is called critical size of flipping, acri- For a in the range 
of 10^^ to 10^^ cm^^, we can estimate Ucri as a function 
of a. 

The obtained results are plotted in Figure [8] when dif- 
ferent internal relaxation processes are accounted for and 
in the presence of stochastic H2 torques. The region 
where grains with a < acri can flip is denoted with flip- 
ping, and the region in which grains with a > acri can 
not flip is denoted with no flipping. 

It can be seen that for a > lO^'^cm^^, acri is similar for 
both paramagnetic and superparamagnetic grains. But 
for a < 10^'^cm^^, acri for the former decreases sub- 
stantially, and reaches very small size ~ 5 x 10~^cm for 
a = lO^^cm"^, while Ucri is about 6 x 10~^cm for the 
later case. Therefore, it exists an area, A, in which the 
superparamagnetic grains can, but paramagnetic grains 
can not flip. 




I 0.00 



.00 
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a/1 0"" cm 

critical size of grain flipping in the ISf' 



Barnett + nuclear + sup relaxation 
Barne:t + n jclear relaxation 
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Flipping for superparamagnetic grainy' 
No Flipping for paramagnetic grains^ 



FLIPPING 



0,0: 



a (cm ^) 

Fig. 8. — Upper: Mean flipping time compared to crossover time 
as a function of grain size for tliree values of active site density, a. 
The grain size at the cross of two timescales defines a critical size 
acri of grain flipping. Lower: Obtained critical size for different a 
for Barnett and nuclear relaxation (dashed line) and total internal 
relaxation (solid line). Regions where grains can and can not flip 
are denoted with flipping and no flipping. 

3.4. Thermal Trapping 

In this section we follow LD99a in describing the prob- 
lem of thermal trapping. This effect happens when the 
flipping is fast, so that the regular pinwheel torques can- 
not appreciably accelerate the grain between the flips. 
As flips happen, the regular pinwheel torques change the 
sign and are averaged to zero. We consider the flipping 
induced by total internal relaxation, in the presence of 
impulses due to II2 formation for the ISM. 

Following LD99a, whether the grain can be spun-up to 
suprathermal rotation after the crossover or not depends 
on the acceleration by the systematic torques and the 
thermal mean flipping time. Let Jtrans be the angular 
momentum at which the acceleration time by systematic 
torques is equal to the mean flipping time. We obtain 



^tf (Jtrans) 



7' 

^ trans 



(26) 



where J^^^ins 



torque given by 



tgas (G- 1)1/2' 

= Jtrans/ Jth, and G is the total systematic 



(27) 
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where ((r^^)^)^/^ is the amphtude of H2 pinwheel 
torques in dimensionless unit given by equation pop . 

When the flipping is rapid, i.e., when thermal trapping 
occurs, the torques in the grain body coordinate system 
change randomly. As a result, the mean angular momen- 
tum is given by LD99a 



J' 



Jth 



th^tf.trap 



(ttf,trap ~\~ ^gas) 



(28) 



where Jtjn as we discussed earlier, is the grain angular 
momentum corresponding to the gas temperature and 
it f, trap is the timescale for the thermal trapping when 
the grain is trapped at a low-J state by frequent thermal 
flipping. From equation (|28p we obtain 



ttf,trap jJ/ JthY ~ 1 
tgas ~ G-{J/Jthf 



(29) 



By setting the mean flipping time tt f obtained from 
calculations for total internal relaxation in the presence 
of H2 impulses (Fig. [7|), equal to that given by equation 
(|29|) . we can obtain solutions for J. Let us denote the 
lower value Ji and the upper one J2. For J < Ji, grains 
are not thermally trapped, and they can be spun-up to 
J = Ji. However, for Ji < J < J2, grains arc ther- 
mally trapped because the thermal mean flipping time 
is smaller than the acceleration time. When J > J2, 
grains are not thermally trapped at all, and able to ro- 
tate suprathermally. Therefore, hereafter, we define J2 
the critical value of thermal trapping. Note that the 
degree of alignment depends on the value of angular mo- 
mentum, i.e on the lower solution Ji, thus, the variation 
of Ji with the pinwheel torques is also important. 

The upper panel of Figure [9] shows ttf/ tgas as a func- 
tion of (J/ Jd)"^ for the ISM and for grain size a = 0.05, 0.1 
and 0.15^m. It can be seen that ttf /tgas tends to in- 
crease with J/Jd] larger grains have ttf /tgas increasing 
faster than small grains with J / Jd- 

In Figure [9] (^ower panel) we show ttf .trap given by (|29p 
for different a and ttf for two grain sizes a — 0.05 and 
0.1/im. The filled circles show the position of J2 and 
Ji. There, it is shown that different grain sizes can be 
trapped at J < J2 with J2 depending on a. Small grain 
a = 0.05/im is trapped at J < J2 — vTTjd even in the 
presence of very strong pinwheel torques (see the lower 
panel of Fig. [H). 

Following LD99a, we can estimate the trapping time 
of grains based on the value of J2 as 



t 



trap 



Jth 



(30) 



Because the trapping size depends on both J2 and a. 
Using equation pO|). we can estimate the trapping size for 
a given trapping time, as a function of a, and the results 
are shown in Figure 1101 The figure shows that grains as 
small as 3 X lO^^cm can be thermally trapped for 2QQtgas 
by the maximal pinwheel torques with a = 10^^ cm^^. 
For the same trapping time, the trapping size increases 
with the decrease of pinwheel torque efficiency. 

Flipping and thermal trapping were invoked in LD99a 
to explain why observations testify that the grains with 
a < 5 X 10^^ cm are not aligned in the diffuse inter- 
stellar gas. Indeed, suprathermally rotating grains are 
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Fig. 9. — Upper panel: Mean flipping time, tff, induced by total 
internal relaxation in the presence of impulses due to H2 formation, 
as a function of (J/Jd)^ for the ISM. The result for gram sizes 
a = 0.05, 0.1 and 0.15/jm are presented. Lower panel: ttf similar to 
the upper panel, and ttf trap ^ ^ function of (J/Jd)'^ for different 
a, and for grain sizes a = 0.05 and a = O.lfim. Filled circles 
depict the angular momenta when the mean flipping time equal the 
acceleration time. For a given a, the decrease of a decreases the 
acceleration time, and thus decreases the thermal trapping range. 



bound to be aligned by the original Purcell (1979) mech- 
anism. The efficiency of the mechanism is somewhat un- 
clear. The original work by Purcell (1979) predicted for 
short-lived pinwheel torques a rather marginal increase 
of the alignment compared with the classical predictions 
in Jones & Spitzer (1967). This conclusion was corrected 
by Lazarian & Draine (1997), who showed that for grains 
larger than a certain size a^ri, high degrees of Purcell- 
type alignment are possible. The size Gcri was estimated 
by Lazarian & Draine (1997) to be around 2 x 10~^ cm for 
the Barnett relaxation. It gets larger than 10~* cm if the 
nuclear relaxation is accounted for (LD99b). Thus, for 
small grains, the original conclusions in Purcell (1979) 
are valid. However, the result in LD99b assumed that 
the thermal flipping time is proportional to the internal 
relaxation time. Using the new result for thermal flip- 
ping, we found that Qcri spams from 10~^ to 10~^ cm, 
which is much smaller than the value obtained in LD99b. 

Thus, if sufficiently strong long-lived pinwheel torques 
dominate the dynamics of small grains, they should ex- 
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Fig. 10. — Trapping size of the ISM grains as a function of 
the magnitude of pinwheel torques for trapping time ttrap = 10, 30 
and lOOtgas . Here we assume that the flipping is produced by total 
internal relaxation in the presence of impulses due to H2 formation. 



hibit degrees of alignment which contradicts to observa- 
tions. The problem dissolves if the torques are either 
weaker than the threshold we discussed above or grains 
are superparamagnetic or pinwheel torques are short- 
lived. In first two cases fast flipping of grains induces 
thermal trapping and grains rotate thermally in spite of 
the presence of suprathermal torques (see LD99a). If 
torques are short-lived and grains are smaller than the 
critical size, the Barnett and nuclear-induced fluctua- 
tions randomize grains during crossovers. In this case, 
the original Spitzer-McGlynn theory should be used and 
it predicts a marginal alignment. 

4. ALIGNMENT BY INTERSTELLAR RATS PLUS 
PINWHEEL TORQUES 

In this section we study the alignment of grains by 
RATs in the diffuse ISM in the presence of II2 torques. 
We consider the influence of long-lived H2 pinwheel 
torques on the RAT alignment when (i) only low-J at- 
tractor points and (ii) both low-J and high-J attractor 
points exist in the trajectory maps. We study for a range 
of grain size from a = 0.02 to 0.2^m. Our calculations 
for thermal flipping showed that the flipping is most ef- 
ficient for the total internal relaxation (i.e., Barnett, nu- 
clear and superparamagnetic inclusions relaxation) in the 
presence of II2 formation (sec §2). Therefore, we adopt 
this flipping rate in our study below. 

4.1. Equations of motion 

The motion of a grain can be separated into the torque- 
free motion of grain axes about angular momentum J 
and the motion of J about the magnetic field B driven 
by external torques, because the former process is much 
faster than the later one (see WD03; HL08a). 

The motion of the angular momentum vector J about 
the magnetic field B subjected to a net torque is com- 
pletely determined by three variables: the angle ^ be- 
tween and B, the precession angle (/) of J around B 
and the value of the angular momentum J (see Fig. [3]). 
The equations of motion of grains by RATs plus II2 pin- 
wheel torques, when the paramagnetic dissipation is dis- 



d(l> -fu^g^^a'^X 
dt 2 Jsin ^ 

^ 7"rad«^^ 

dt 



2J 



(G(C,0,^,J))-f^B, 
(F(C,0,^,J)), 



dJ 1 



J 



(31) 
(32) 
93) 



-17 = ^l^Tad"- ^(^(^' ■^)) + -cos 

uX Z '-'gas 

where (i;^(f, V, J)), (G(C, 0, ^, J)), and (iJ(C, 0, ^, J)) 
are the RATs components averaged over the torque-free 
motion and thermal fluctuations (see HL08a); ^Ib is the 
Larmor precession rate, a is the grain size, 7 is the degree 
of anisotropy of radiation, A is the mean wavelength of 
radiation field, Urad is the radiation energy density, and 
tgas is the gas damping time (see Table 1). Here cos 9 oc 
J cos 9f{6, J)sm 6d6 with 6 angle between ai and J and 
/(0, J) is the thermal distribution function. is the 
component along ai axis of H2 pinwheel torques. 

For the ISM, the Larmor precession rate is larger than 
the gas damping and the alignment rates; and thus, we 
can average equations p2|) - ((33)) over a precession period 
for (f) from to 27r. As a result, equations ([3T|) -([33 | can 
be reduced to a set of equations for ^ and J, in which the 
spinning and aligning torques (F) and (H) are replaced 
by (F)^ and {H)^, resulted from averaging correspond- 
ing RATs components over the Larmor precession angle 
(p}^ In dimensionless units (J' = J/Jthit' = t/tgas) and 
taking into account the grain flipping, we obtain 



d^ M 



-((F)^,+/+ + (F>^,_/_), 



dt' J 

dJ' 

— ^M{{H)^,+f+ + (i/)^,-/-) - J' + r'^,cos eif+ 

where is the dimensionless H2 torques. 



(34) 



'i.Jth ' 



(36) 



and /_(_ , /_ are the fraction of the time interval dt the 
grain stays in the positive and negative flipping states, 
respectively (see Appendix D); (F)0(+ _) and (i?)0(+__) 
correspond to positive and negative flipping states (see 
HLOSa).!-* 

We adopt the RAT efficiency Qr from the AMO with 
the ratio of two first components of torques q™-"-^ = 1.2 
(see Appendix C). We first average RATs over torque- 
free motion, and thermal fluctuations for the irregular 
grain with ratio of inertia moments Ii : I2 '■ I3 — 1.745 : 
1.61 : 0.876 (see WD03; HL08a) to get (Qe^). Then, we 
calculate aligning and spin-up components of RATs us- 
ing equations (C13) and (C15) where Qe ,j = 1 — 3 are 
replaced by (Qcj)- The resulting aligning and spin- up 
components of RATs after averaging over precession an- 
gle ip are shown in Figure [Tl] for three values of J and 
ip = 0°. For this setting of AMO, the two stationary 

For the sake of simplicity, hereafter, we denote (F) = 
{F{i,<l>,i,,J)), (H) = {H{i,<t>,^l,,J)), and {F)^ = 0, ^, J)>^, 

The positive and negative flipping states are defined based on 
the initial directions of ai with respect to J. The positive flipping 
state corresponds to ai paraUel to J and the negative flipping state 
corresponds to ai antiparallel to J. 
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Fig. 11. — Aligning (upper) and spin-up (lower) components 
of RATs for the ellipsoidal AMO witli g™"^ = 1.2 and ratio of 
inertia moments Ii ■ I2 '■ I3 = 1.745 : 1.61 : 0.876 and for light 
direction 1/1 = 0°, obtained by averaging over thermal fluctuations 
for different values of angular momentum J. Lower J corresponds 
to stronger thermal fluctuations and the magnitude of RATs de- 
crease. Functional forms of RATs also change as J decreasing, but 
stationary points ^ = and n remain unchanged. 

points COS = ±1 are expected (see the upper panel). 
For the spin- up component, (H)^ at these stationary 
points changes the sign to the opposite when J ^ 1. As 
a result, we expect low-J attractor points in trajectory 
maps for this alignment (see HLOSa). 

4.2. Trajectory maps 

Using RATs from Figure [11] we solve equations of mo- 
tion for J and ^. We use initial conditions Jq = Jth and 
^0 is generated from a uniform distribution of the angle 
between J and B. For long-lived H2 torques, we assume 
t']^ = 10^ and a constant time step St' — 10~^. We as- 
sume again that H2 torques are parallel to the axis of 
maximal inertia for positive flipping state. 

The trajectory maps with RATs from Figure [TT] are 
presented in Figure 1121 The upper panel represents the 
alignment by RATs only, and the lower panel repre- 
sents the alignment by both RATs and the long-lived 
H2 torques. 

The upper panel of Figure [T^] shows the trajectory 
maps for the alignment by RATs only. For this case, 
all grains with different initial angles ^ are driven to a 
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Fig. 12. — Phase trajectory maps for the alignment in the diffuse 
ISM of a O.l^m grain for light direction ip = 0°i in the absence 
of H2 torques (upper panel) and in the presence of long-lived H2 
torques with a = 5 X 10^^ cm~^ (lower panel). RATs align grains 
at a low-J attractor point B seen in both panels, but the presence 
of long-lived H2 torques gives rise to the new high attractor point 
(A) corresponding to suprathermal rotation in the lower panel. The 
AMO with torques ratio q"^'^^ = 1.2 and the ISRF are adopted. 

low-J attractor point A. In the presence of H2 pinwhecl 
torques, the lower panel of Figure [12] shows that such 
long-lived torques can produce a new attractor point of 
high angular momentum (point B) with J/ Jth ~ 9. The 
alignment angle corresponds to the perfect alignment of 
J with respect to B (i.e. cos f = —1). 

The appearance of a new high-J attractor point at 
cos ^ = — 1 in the presence of long-lived H2 torques can 
be easily understood. First, some grains with initial di- 
rection in the vicinity of cos ^ = are spun up by pos- 
itive spin-up torques (see Fig. Ilip so that flipping rate 
decreases. Thus, pinwheel torques become more efhcient 
and act to increase J. When J > 3 Jth, the flipping 
decreases substantially and the grains are able to escape 
from the thermal trapping (grains with trajectories 1 and 
2). At the same time, the aligning component of RATs 
drives the grains to the stationary point cos = — 1 as 
in the absence of pinwheel torques. Finally, these grains 
are stably aligned at the high— J attractor point. The 
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rest of grains are aligned at low- J points. 

In fact, as seen in Figure [HI the aligning torque has 
two stationary points at cos = ±1 for J Jth- 
In the absence of II2 torques, this stationary point is 
a low-J attractor point because of the spin-up compo- 
nent of RATs {H)^{^s = tt) < 0. In the presence 
of positive H2 pinwheel torques, which are assumed to 
be larger than RATs, we have the total spin-up torque 



{H)tot = 
because 



M{H) 



rj > 0. Thus, (H) 



< 



< (see Fig. [TT])- This indicates 

-TT 

that the stationary point cos i^s = — 1 becomes an attrac- 
tor point. 

We stress that the role of II2 torques in producing the 
perfect alignment of J with B in Figure [T^ docs not in- 
volve the paramagnetic dissipation. Instead, the station- 
ary points are produced by the RATs (see LH07a). 

4.3. Increase of J at low attractor points 

In the absence of pinwheel torques, the value of J at 
low-J attractor points (hereafter Jiow) is constrained by 
the internal relaxation as well as the negative spin-up 
component of RATs. The presence of pinwheel torques 
increases the value J/^^ as (see eq. inS|) 



= -ej,e^j'Zi^ + jZl^{f+ - /-), (37) 

where we used the fact that ^ {H)^- and /+ + 

/_ = 1 in eliminating the term /+ , /_ associated with 
(iJ)^, ej/, evaluated at J' — J'l^w^ ^^'^ factor by which 
{H)^ is reduced by averaging over thermal fluctuations 
at fixed J', is a function of ip describing the decrease of 
spin-up torque with ip. At J — Jth, we have ej ^ 0.1 (see 
HL08a). From our calculation of RATs using DDSCAT 
for the shape 1, we found that = 1 for ij) = Q° and 
~ 0.25 for = 70° for a 0.2^m grain (see also HL08b). 

From equation (|57|) it is obvious that for fast flipping, 
i.e., /+ — /_ = 0, the pinwheel torques play a marginal 
role in spinning up grains. However, our calculations in 
§2 show that the flipping gets slower for grains larger 
than 10~^cTO, so the pinwheel torques can have an im- 
portant effect for these grains. 

Using equation (jD4p and adopting the time step dt = 
IQ^^tgas, we can estimate /+ — /_ for different J and 
grain sizes assuming the grain is in positive flipping state 
at the beginning of the time interval. Figure [13] shows 
that for grain smaller than ^ 3 x lO^^cm, we have /+ — 
/_ « 0. For larger grain a — lO^^cm, /+ — f- — 0.04 for 
J = Jth = i-SJd- Combining equations (fT0|) . (fTSjl and 
(|37|) with Figure [TS] we can obtain J/^^ as a function of 
a and grain size. 

We note that the mean flipping time increases with J, 
therefore, at some value of J (i.e., J > 5), the grain does 
not flip. Thus, /+ — /_ = 1, and the low-J attractor 
point transfers to high- J attractor point (see §4.2). 

To get the exact value of J'l^^, we solve the equations 
of motion and construct the trajectory maps. Figure 
[T^ represents the obtained results using the AMO with 
torques ratio g""*^ = 1.2 for two radiation directions 
V' = 0° (upper panel) and 70° (lower panel) and for three 
values of grain size (a = 2 x 10~^, 6 x 10^^ and 10^^ cm). 
For both directions, the value Jiow increases fast with 
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Fig. 13. — Difference between the fraction of the time the grain 
spends in the positive flipping state and the fraction of time the 
grain spends in the negative flipping state as a function of the grain 
size for J = l.SJd and 2.5Jd assuming that the grain is in positive 
flipping state at the beginning of the time interval with time step 



dt = io~-^t„ 



f+- f- < 10"'" for a < 0.3 X 10" 



the decrease of a, i.e., the value of II2 torques increases. 
However, for the grain with a < 2 x 10~^ cm, Jiow can 
not exceed the thermal value even with the strongest 
H2 torques. This arises from the fast flipping of such 
small grains that decreases significantly the effect of H2 
torques. For the former case, Jiow < Jth for all possible 
magnitudes of pinwheel torques. This result is consistent 
with our calculation in the preceding section. 

The above results are obtained using the default el- 
lipsoidal AMO. However, since the H2 torques depend 
weakly upon the shape of grains, we expect the results 
are valid for irregular grains. 

4.4. Effect of pinwheel torques on the alignment with 

high-J 

Now let us consider the effect of pinwheel torques 
on the RAT alignment in the ISM for a large grain 
a — 0.2/im for which the flipping is inefhcient, and the 
torque ratio = 0.78 to get the alignment with high- 
J attractor point for the radiation direction tp = 70° (see 
LH07a). We assume that the pinwheel torques are paral- 
lel to ai axis and have a magnitude with a = lO-'^^cm"^. 

Figure [T^] shows the trajectory maps in the absence 
and presence of pinwheel torques. In the first case, RATs 
align grains with a high-J attractor point A and a low- 
J attractor point B. In the second case, the presence of 
pinwheel torques parallel to ai axis increases the angular 
momentum of the attractor points A and B. The low-J 
attractor point becomes the high-J attractor point, but 
the alignment angle cos ^ = —0.8 does not change. 

When the amplitude of pinwheel torques is less than 
that of RATs, the pinwheel torques can provide a new 
type of attractor point. Indeed, at a high attractor point 
the pinwheel torques may act opposite to RATs, mov- 
ing point A to the position D. The point C is a regu- 
lar crossover point with the crossovers assisted both by 
RATs and pinwheel torques. As the crossover happens, 
the grain flips (see Spitzer & McGlynn 1979, Lazarian 
& Draine 1997, 1999a) and the RATs and the pinwheel 
torques start acting in the opposite directions. This cre- 
ates a new uplifted low attractor point similar to B, but 
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Fig. 14. — Value of angular momentum at low attractor points 
as a function of the active site density a for different grain sizes 
and two light directions 1/1 = 0° {upper panel) and ip = 70° {lower 
panel). The AMO similar to Fig. I12l is adopted. 

corresponding to the opposite state of rotation. The in- 
dividual realizations of pinwheel torques make the trajec- 
tory maps not symmetric, but if we consider an ensem- 
ble of grains with equal number of positive and negative 
pinwheel torques the trajectories for AMO get symmet- 
ric and we can show only an upper part of them as we 
do in Figure [TSl 

The most important effect of the strong pinwheel 
torques is their possible uplifting of the low attractor 
point. This should increase grain alignment in the case 
when the trajectory maps do not have high attractor 
points. The latter situation corresponds to a substantial 
portion of the parameter space (see LH07a). 

4.5. Role of RATs in spinning-up grains 

In §5.2 we found that positive RATs can support the 
pinwheel torques to bring grains to high- J attractor 
points because they are not fixed in the grain body axes. 
Now let us estimate the value of J that spin-up RATs 
can induce the grain for a given light direction ^. 

The spin-up component of RATs for a light direction 
^ = 70° is shown FigureE] For J > Jth > Jd, {H)^ > 
for cos ^ > 0, and (H)^ < otherwise (see also the lower 
panel of Figlll|). The magnitude of spin- up torque at the 
angle ^ = is given by equation (|12p when the radiation 
is parallel to the magnetic field. For an arbitrary radia- 
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Fig. 15. — Trajectory maps of the RAT alignment for a grain size 
a = 0.2iJ,m without and with pinwheel torques (upper and lower 
panels) for a = 10^^ cm~^. The presence of pinwheel torques par- 
allel to ai axis results in the high-J attractor point B, but does not 
change the alignment angle. C and D are expected attractor points 
when pinwheel torques are anti-parallel to ai and the amplitude 
of pinwheel torques is smaller than that of RATs. 

tion direction ip, the magnitude is decreased by a factor 
of . We can estimate the maximal angular momentum 
spun-up by RATs when radiation makes an angle ip with 
B as foUowings 
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where J^ax^ is given by equation 
For A > 1.8a, rj — —3, we obtain 
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It can be seen that J, 



(39) 



size as a function of a 
obtain J' „, 



increases rapidly with the grain 
For a grain size a_5 = 2, we can 
12 for ip = 70° with e.0 = 0.25 and e./^ = 



7/2 



1. When grains are spun-up by RATs to Jmax '3> Jth, 

Wc note that as J' increases, the factor ej/ in equation Il39p . 
which describes the decreases of the angular momentum due to the 
thermal fluctuations, increases, and achieve the order of unity for 
J' > 10 (see HL08a) 
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flipping stops, and H2 torques become efiicient. From 
equation ([39]) it can be seen that for grains a-5 > 1, 
Jmax — 1 every ip. This indicates that RATs for the 
ISRF are able to spin up large grains to suprathermal 
rotation, and thus they are not thermally trapped. 

5. DISCUSSION AND SUMMARY 

5.1. Grains are flipping and get trapped 

New elements of grain dynamics, namely, thermal flip- 
ping and thermal trapping were reported in LD99a. The 
validity of these findings was challenged in Weingartner 
(2008), who found that grains can not experience thermal 
flipping as a result of internal relaxation. In the present 
work we showed that although internal relaxation can 
not result in grain thermal flipping, but in the presence 
of impulses due to H2 formation and gas bombardment, 
grains can flip, and get trapped. The impulses play the 
role of an engine to bring grains through the boundary 
ai _L J. 

The thermal trapping is an essential process that, ac- 
cording to LD99a, explains why, in accordance with ob- 
servations (see Kim & Martin 1995), small grains are not 
aligned by the Purcell (1979) paramagnetic relaxation 
process. The alternative explanations to this observa- 
tional fact (see Lazarian 1995) are much less appealing. 
The sizes of the thermally trapped grains are not differ- 
ent from those which can be obtained using the LD99a 
treatment, but smaller than those discussed in LD99b. 
As a result, pinwheel torques are more important for in- 
terstellar grains than one could infer from LD99b. 

5.2. Extending theory of crossovers 

Crossovers are times when grains subject to pinwheel 
torques slow down, then flip and get accelerating. The 
original theory of crossovers was suggested by Spitzer & 
McGlynn (1979). Lazarian & Draine (1997) improved 
the theory by accounting for the value of the resid- 
ual thermal angular momentum, associated with ther- 
mal wobbling of the grain (see Lazarian 1994). LD99b 
showed that with nuclear relaxation taken into account 
the theory for regular crossovers in the spirit of the afore- 
mentioned works is applicable only for grains larger than 
~ 10"** cm, which excludes most of grains in diffuse in- 
terstellar gas. LD99a, instead, introduced the concepts 
of thermal flipping and trapping for smaller grains. 

The present paper shows that the flipping and trapped 
grains are smaller than those discussed in earlier papers. 
Therefore, in the absence of RATs, there is a range of 
grains for which the crossovers happen over the time 
scale larger than the internal relaxation time. In this sit- 
uation the thermal value of the angular momentum can 
decrease below a result of the action of the pinwheel 
torques. This increases the randomization of grains dur- 
ing crossovers for the range of sizes larger than the flip- 
ping size. 

5.3. How common is suprathermal rotation? 

Suprathermal, i.e. much faster than thermal, rotation 
was assumed to be default for most of the interstellar 
grains after the Purcell (1979) study. Studies of RATs 
in DW96, DW97 only made this point stronger. How- 
ever, later works that reported thermal trapping of grains 
(LD99a, LD99b) and low-J attractor points (Weingart- 
ner & Draine 2003, LH07a) made one wonder whether 



most grains rotate thermally. The study in HLOSa re- 
vealed that in the presence of gaseous bombardment or 
other randomization processes, the grains tend to diffuse 
from the low-J to high- J attractor point, provided that 
such high-J attractor point coexist with the low-J one. 
According to the parameter space study in LII07a (see 
Fig. 24 in LII07a) this corresponded to a fraction of 
grains that could vary depending on the angle between 
the radiation direction and the magnetic field, as well as 
on the g™"^ ratio. 

In the present study we showed that for sufficiently 
large grains the pinwheel torques, provided that they are 
comparable with or stronger than the RATs, can induce 
suprathermal rotation of grains. Interestingly enough, 
this can happen even in the situations when only low-J 
attractor points exist. Therefore, the presence of pin- 
wheel torques can increase the percentage of the grains 
rotating suprathermally. However, it is worth noting, 
that for a substantial portion of the parameter space ra- 
diative torques act against pinwheel torques attempting 
to decrease the grain rotation rate. 

5.4. Implications to polarization modeling and magnetic 
field diagnostics 

A quantitative theory of grain alignment is very im- 
portant to polarization simulations. A number of works 
dealing with polarization simulations assume the perfect 
alignment of the grain axis with respect to the magnetic 
field (see CL07; Bethell et al. 2007; Pelkonen et al. 2007; 
Falceta-Goncalves, Lazarian & Kowal 2008). A detailed 
comparison between the polarization efficiency predicted 
by the theory of grain alignment and observations of 
dense cloud in Whittet et al. (2008) shows that a high 
degree of alignment is required to explain the observed 
polarization degree. How can this be accomplished? 

In Lazarian & Hoang (2008) we proposed one way of 
increasing the degree of grain alignment, namely, we no- 
ticed that in the situation when grains have superparam- 
agnetic inclusions the high-J attractor points are present 
for all the angles between the radiation direction and the 
magnetic field, as well as on the ratio. Thus the 

gaseous bombardment is expected to move all grains to 
the high attractor points, making the alignment perfect. 
In the present paper we show an alternative way of in- 
creasing the degree of grain alignment. In particular, 
we show that RATs plus pinwheel torques can enhance 
the degree of alignment by increasing the value of angu- 
lar momentum at the low-J attractor point. However, 
the circumstances for which the pinwheel torques affect 
the RAT alignment are rather restrictive: the pinwheel 
torques are required to be comparable or stronger than 
RATs. 

With these theoretical results one can attempt to dis- 
tinguish between the two alternatives above and simul- 
taneously get insight into the grain composition. First 
of all, more studies similar to those done by Whittet et 
al. (2008) arc required to make sure that the high align- 
ment is not a result of favorable illumination geometry. 
Alternatively, one may make measurements in situations 
where the illumination geometry is known. Comparing 
these results with the analytical and numerical predic- 
tions (see LH07a, HLOSa) one can establish with higher 
certainty that the grain alignment is enhanced compared 
to what is expected from ordinary paramagnetic grains 
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subject to RATs in the absence of the pinwheel torques. 
In paraUel, one can search for the correlation of grain 
alignment enhancement with the environments where we 
expect the higher amplitudes of the pinwheel torques. 
If such an enhancement is present, it will indicate that 
grains are not superparamagnetic. If they were super- 
paramagnetic, we would expect the perfect alignment for 
typical conditions of the diffuse ISM, which one cannot be 
improved further by the action of the pinwheel torques. 
No correlation, but high degrees of alignment mean, on 
the contrary, the existence of superparamagnetic grains. 

One can conclude by observing an interesting cycle in 
the development of the grain alignment theory. The ini- 
tial theory by Davis-Greenstein (1951) appealed to para- 
magnetic dissipation in ordinary paramagnetic grains. 
The mechanism was shown not to be sufficiently strong, 
however. Then to enhance the efficiency of the ini- 
tial Davis-Greenstein process both superparamagnetism 
(Jones & Spitzer 1967) and pinwheel torques (Purcell 
1979) were appealed to. At the moment we may face 
the situation that the RATs are not strong enough to 
explain the alignment and again we appeal to super- 
paramagnetism (Lazarian & Hoang 2008) and pinwheel 
torques (this paper). Nevertheless, there is a substantial 
difference between the two situations. The paramagnetic 
relaxation in the diffuse interstellar medium could pro- 
vide a few percent alignment at most. The RAT align- 
ment provides the degrees of alignment about 20% even 
in the cases of unfavorable directions of illumination (see 
HLOSa). In addition, while the Davis-Greenstein process 
favored the alignment of small grains, the RAT align- 
ment, in accordance with observations, favors the align- 
ment of large grains. 

5.5. Summary 
Our main results are summarized as follows: 

• We extended the definition of the pinwheel torques 
to include three new processes, namely, radia- 
tive torques arising from grain infrared emission, 
torques due to interactions of grain with electrons 

in plasma, as well as, torques arising from gas flow 
interacting with a helical grain. We proposed to 
consider the radiative torques arising from the ac- 
tion of isotroptic radiation on a grain as pinwheel 
torques. All these torques we identified as long- 
lived, as their life-time is expected to be longer 
than the typical grain rotational damping time. 
We used II2 torques as a proxy of for general pin- 
wheel torques, but, unlike earlier studies, we al- 



lowed the pinwheel torques to be both short-lived 
and long-lived. We found that in the ISM, the pin- 
wheel torques due to infrared emission is in general 
weaker than that due to H2 formation, and RATs. 

• We proved that flipping and trapping of astrophys- 
ical dust grains is the active and important pro- 
cess. We did this by augmenting the torques arising 
from internal relaxation by the stochastic torques 

that inevitably accompany the action of pinwheel 
torques and gas bombardment. . We calculated 
the thermal flipping induced by the Barnett, nu- 
clear relaxations, as well, as internal relaxation in 
superparamagnetic grains in the presence of im- 
pulses due to H2 formation and gas bombardment. 
Adopting the revised diffusion coefficients for inter- 
nal relaxation from Weingartner (2008), we foimd 
that flipping is fast for small grains, and increases 
when superparamagnetic inclusions are accounted 
for. We obtained the critical size of flipping and 
trapping size of the thermally trapped grains for 
the ISM. 

• The most important result of this study is related 

to the increase of the expected degree of grain align- 
ment when dust grains are subject to both RATs 
and pinwheel torques. This increase stems from the 
increase by pinwheel torques of the value of the an- 
gular momentum at the low- J attractor point. The 
increase of angular momentum decreases the wob- 
bling of grains at the low- J attractor point and 
therefore increase the degree of internal alignment 
as the internal alignment (of ai with J) is nearly 
perfect when the grain rotates with rates much 
larger than the thermal rotation rate. The joint 
pinwheel plus RATs alignment is different from the 
Purcell (1979) alignment. In the latter mechanism 
the alignment is due to paramagnetic dissipation, 
which induces torques that are very weak compared 
to RATs. 
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APPENDIX 
DIFFUSION COEFFICIENTS 

LR97 derived the diffusion coefficients for the Barnett effect by assuming that in thermal equilibrium, the probability 
current 

S = A{e)fTE{0,Td) - ~ {B{6)fTE{e,m . (Al) 

with fTE{{S,Td) the Boltzman distribution function, is equal to zero. 
Weingartner (2008) used instead of angle 9 a variable q = ^^j^. where for an oblate spheroid E = 



and found diffusion coefficients different from those in LR97, namely 
(q-m-q) 



A{q): 



tint (J) 



(A2) 
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+c{c)ih - qy/^cMCq) - c^'^ [3 + 2C(/i - 1)] 

x{h- g)^'^exp[-C(/i - q)] / exp(x^)(ix, 



(A3) 

where h — /1//2, C — (J/Jd)'^, B{q — h) — 0. We would claim that the main difference between these coefRcients and 
those in LR97 is that both A and B coefficient get zero when the angle is equal 9 is equal to 90 degrees. 

Dealing with oblate grains, we prefer, however, to study the dynamics in terms of 6. The corresponding Ito's formulae 
(see Gardiner 1983) for the transformation of variables in the stochastic calculus provide us with: 

de . [A[q)- + + v/5(^^rf-, (A4) 

where dw is a random variable with {duP') — dt. Thus, using the equivalence of Langevin equation and Fokker-Planck 
equations, we obtain 

^rn^ ^r^^e B{q) 8^0 -{h - 1) . ^ ^ , -COS 26* 

oq 2 oq'' 2Unt(J) (/i — l)"^sm 20 

B{e) = B{q)(^^y , (A6) 

where we used the relation q — 1 + {h ~ l)sin^0 to calculate first and second derivative of 9 with q. Equation (|A5I 
differs with what in Purcell (1979) and LR97 by an additional term arising from variable transformation. 

However, while considering grain dynamics, we introduce, in addition to the regular external torques, stochastic 
external torques. These torques are inseparable part of the grain dynamics. Following Lazarian & Draine (1997), 
we consider both stochastic torques arising from H2 formation and those arising from gaseous bombardment. Our 
treatment can be easily generalized for a different nature of the pinwheel torques. The diffusion coefficients due to all 
active sites on the grain surface in the grain-body system i, y and z, are given by 

({AJ^f) ^-fHUHa^vnaVa^^^, (A7) 

((AJ,)2) = ((AJ.)2), (A8) 

{(AJ,)^) ^jHriHa^vna^plr^, (A9) 

where 7/7 is the fraction of H atom converted to H2, nn is the gas density, vh is the velocity of incoming atom, 
and Tgas is the gas temperature, and r||,rj_ are geometrical factors depending on the grain shape (see RDF93). For 
simplicity, we assume Fji = r_L = 1. In dimensionless units, 

((A4)2) = (AlO) 

{{AJ'^f) = ((AJ;)^), (All) 

{{AJ',f) = -^^, (A12) 

where Ekin — Po/'^'^nH is the kinetic energy of evaporating H2 molecule, and we neglected the anisotropy of diffusion 
coefficients. Here we adopt Ekin = 0.2eV as suggested by Purcell (1979). Transforming the diagonal tensor B'^j with 
B'^^ = ((AJj'j)^), and Bl_^j = from the grain body system to the spherical system described by J, 9, (j>, we obtain the 
component along 9 axis 

BH2i0) = 0.5(5^1 + 5^2) (sin'*6l + 0.5cos26' + sin26') 

+1.5B33(sin26lcos26l), (A13) 

with = = ((AJ;)^), 5^3 = ((AJD^)- 

The fluctuations arising from gas bombardment contribute also to the grain flipping. The diffusion coefficients 
obtained in Roberge et al. (1993), in the zeroth approximation, are given in the grain body system of reference by 

1 / T 



BIT = -nnmvHkTgasb'' [l + {n - 1)^ j (A^ + F^), (A14) 
Bfr-S^f, (A15) 
Bit = \nHmvHkTgasb^ + (« - 1)^) (A|| + F,,), (AI6) 
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where uh^vh are density and velocity of atomic gas, respectively, Tgas,Td are temperature of gas and dust, k is the 
fraction of H atoms sticking to the grain that evaporate, A|| and are geometrical factors of order of unity. 

In dimensionless units 5^™" = -^if^ Jth/^gasi we obtain 



" = (^1 + ('^ - 1)^ j (Ax + Tj.), (A17) 
52^2°" -Si™", (A18) 
= (l + - 1)^) (A|| + r,,) (A19) 

FLIPPING FORMALISM 

Below we present the Fokker-Planck equation approach to study the flipping dynamics of grains. We note that this 
approach is presented in the classic book of Gardiner (1983), and brought to grain dynamics context by Roberge & 
Ford (2000; RFOO).!^ 

Consider a brick grain with sides b x b x a and a < b. The energy of this grain with temperature Td is given by 

m = ^ (1 + - ' (Bi) 

where Ii is the moment of inertia along the maximal inertia axis ai, Td is the dust temperature, 9 is the angle between 
ai and the angular momentum J of magnitude J, ks is the Boltzmann constant. Thermal fluctuations induce changes 
of and when 9 > n/2, the grain undergoes a flipping. 

We assume that at time t — 0, the grain is present at a given position 9 £ [0, 7r/2]. We can ask, after a time interval 
t, what is the probability that the grain is still present in the initial range. We denote such a probability by G{9,t). 
After the time t, the grain can be present at position 9' G [0,7r/2] with the probability f{9',t), so that we have 

G{9,t)= r'%{9',t)d9'. (B2) 
Jo 

We use the conditional probability p(02) i2|^ij i-e., the probability that the grain present at 9i,ti will be present 
at 6*2 at the time ^2- In terms of p(02, ^216*1, ii) 

f{9',t)^p{9',t\9,Q)^p{9'M9,^t), (B3) 

where the conditional probability which is uniform in time is adopted. 
Hence, the probability of flipping from the range [0, 7r/2] to [7r/2,7r] is 

To find G{9, t) from equation (|B2p . we need to know f{9', t). The initial condition t = states that 

f[9',t = Q) = 5{9' -9). (B5) 

Thus, 

.7r/2 i.tt/2 

G{9,t)^ f{9',t)d9'= p{9',0\9,~t)d9'. (B6) 
Jo Jo 

It is well known that if / satisfies the Fokker-Planck equation, then the probability p with initial condition satisfies 
the backward Fokker-Planck equation: 

dp{9',0\9,t) dpi9',0\e,t) Beg(Pp{(mm 

dl "^'^ d9 2 W • ^^^^ 

Integrating equation (|B7P with 9' from to 7r/2 combined with equation (|B6|) and then changing variabel t — > —t, we 
obtain 

dG{9,t) ^ dG{9,t) , Bge d^G{9,t) 

-^ = "^'^9-^-2 — m^- (^^^ 

Initial and boundary conditions for equation (|B8P are given by 

G{9, i = 0) = 1 for < 61 < 7r/2, (B9) 
dG 



dt 



= 0. (BIO) 

t=o 



^® We provide more details of the RFOO approach, as this paper has circulated only in the form of the preprint and may not be readily 
available. Their treatment is more rigorous than that in LD99a. 
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Also, a condition for the grain to flip immediately at 6* = 7r/2 is 

G(6i = 7r/2,i) = 0. (Bll) 

This boundary condition corresponds to an absorbing boundary, i.e., the grain "escapes" at 9 = 7r/2. The problem 
of a particle escape with initial and absorbing boundary conditions (jBlip is discussed in a great detail in Gardiner 
(2004; §5.2.7). 

The mean mean flipping time is defined by means of the flipping probability (eq. |B4|: 



ttf{0)= f t<^{t)dt^~-( t^^^^^dt ^ -tG\^ + f G{9,t)dt^ f G{e,t)dt, (Bf2) 
Jo Jq Jo Jo 

where we used G{6,t — oo) = 0. By integrating equation IB8I over time and taking use of equation (|B12p . we obtain 



Hence 



d{Uf){9) Beg dHt,fm 

TTn + ^ ■ (^14) 



86 2 86 

The initial and boundary conditions for equation (|Bf 4p read 

Uf{e = ^/2) = 0, (B15) 

^(^? = 0) = 0. (B16) 

The boundary condition (|B16p implies that the grain flips immediately when the maximal inertia axis perpendicular 
to the angular momentum (i.e., 7r/2 is an absorbing boundary). 

The solution of equation (|B14p can be given by (Gardiner 1994; RFOO) 

where 

v&(0')=exp(^£|^d0"^. (B18) 

AMO AND RATS COMPONENTS 

Here we present RATs produced by a radiation beam on a helical grain as shown in Figure [2l and the averaging of 
RATs over thermal fluctuations. 

CI. RATs and mean values over the spectrum of radiation fields 
The radiative torque resulting from the interaction of radiation field with a grain of size a is defined by 

i-rad — ^ *4rJ \'^^) 

where 7, A, and Urad are the degree of anisotropy, mean wavelength and total energy density of radiation field, 
respectively, Qr is the radiative torque efficiency vector, and overlines denote the averaging over the spectrum of the 
incident radiation field. 

Let ux be the energy density per unit wavelength A, the mean wavelength and total energy density for the radiation 
field is defined as 

<-) 

Urad = J UxdX. (C3) 

Let Qr be the radiative torque efficiency produced by the monochromatic radiation field with wavelength A, the mean 
torque efficiency over the radiation field is 

Q-r = l^^. (C4) 
J u\d\ 
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C2. General expressions of RATs 
RAT for the toy model in Figure is given by the similar form as equation (jCip : 

r™d=2^^^^Qr, (C5) 

where I2, defined as the grain size a, is the size of the squared mirror. 

Using the self-similar scaling of the magnitude of RATs obtained for irregular grains of size a induced by the radiation 
field of wavelength A, 



|Qp|^0.4[^-j for A > 1.8a, (C6) 

-0.4 for A < 1.8a, (C7) 
and the functional forms of RATs from the AMO, we can write RAT components as following 

O IQrk— ge3(9,/3,'i> = 0) . . 



(7e,(e, /?, $ = 0) = Lma^^^^^^^J + ^cos /?sin 29 - ^cos /3sin 29 - ^cos 2fi\ , (Cll) 



where 



(7e2 (9, /?, $ = 0) = — !-C cos /3cos ^9 ^^cos ^/3sin 29 ^-^sin 29 + n^cos /3sin ^9 , (C12) 

(7e3 (9, ^, $ = 0) = ^Cnisin /3 [nicos 9 - nzcos /3sin 9] + f ^) ^ (s^ _ l)X(9)sin 29, (C13) 
A V'2/ A 

with = max {Q , (3 , ^ = for j = 1,2 and 3. The magnitude ratio of torque components g™"^ is defined 

by 

^ax ^ mag;(Qei(9,/3,$ = 0))/3 

maa:(Oe.(9,/3,$ = 0));3' ^ ^ 

In equations (|Cll[l - (jC13p . C is a function given by 

C = |?iiCos 9 — n2sin 9cos /3| , (CIS) 

where 9 is the angle between the axis of major inertia ai and the radiation direction k, /3 is the angle describing 
the rotation of the grain about ai (see Fig. ^ower); ni = —sin i, 712 = cos i are components of the normal vector of 
the mirror tilted by an angle i in the grain coordinate system, a, b are minor and major semi-axes of the spheroid, 
s = a/b < 1 and e is the eccentricity of the spheroid, li is the distance from the mirror to the spheroid, and I2 is the 
size of the squared mirror; K{Q) is the fitting function (see also LH07a). The second term of equation (jC13p represents 
the torque due to the spheroid. Assuming li ^ X and b, a ~ ^2 ^ ^ii then this term is subdominant compared to the 
first term. Thus, we disregard it in our calculations (see Fig. 

Our calculations for the alignment in the presence of thermal fluctuations showed that the AMO can reproduce the 
alignment property with low-J as found with RATs obtained by DDSCAT when q^^ is modified to (see HL08a) 

4/1 / f3rOS ^9 — ll 77^ 77^ \ 

qe, (9, /5, $ 0) -iC ni7i2^ + ^cos /3sin 29 - ^cos /3sin 29 - 77in2Cos 2/3 . (C16) 

A V 2 2 2 / 

This modification can arise from the imperfect scattering and/or the absorption effect by the mirror (LH07a). Also, 
the results in LH07a remain unchanged because the averaging over f3 for the last term goes to zero. 

We adopt the AMO with i = 45° in this paper, unless mentioned otherwise. We also assume the amplitude of Qe^ 
is comparable to that of Qei and • 

RATs at a precession angle $ (see Fig. ^ower) can be derived from RATs at $ = using the coordinate system 
transformation, as follows: 

Qei(e,/3,$)-ge,(9,/3,<i> = 0), (C17) 
Qe, (9, /?, $) = Qe, (9, /?, $ = O)cos$ + (6, /3, $ = 0)sin$, (C18) 
Qe, (9, /3, = ge2(9, /3, $ = 0)sin$ - Q^, (9, /3, $ = O)cos$. (C19) 
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To study the alignment of the angular momentum with respect to magnetic field, we use the spherical coordinate 
J, ^ and (f) (see Fig. [3|). In this coordinate system, RATs components are given by 

F{tp, (/), —Qei {£., V'l sin tAcos fcos <j) ~ sin ^cos -0) + Qeii^, "0; 0)(cos ipcos ^cos (/) — sin ^sin ^p) 

+Qe3 , V', 0)cos ^sin </., (C20) 
G'(V', (/>, C) = Qei , V', '/')sin V'sin </> - Qea (Ci V', '/')cos -0sin + '!/'> «^')cos 0, (C21) 

H{ip, (/), ^) =Qei {£., V'l ?!')(cos tAcos ^ — sin ipsin ^cos 0) + Qe-^iS,, tp, 0)(sin'0cos ^ + cos V'sin ^cos 0) 

+ge3(e,i/',0)sin ^sin 0, (C22) 

where Qei (^, "07 Qe2 (^j "01 0)j Qea (Ci 01 0): ^-s functions of ^, -0 and 0, are components of the RAT efficiency vector in 
the lab coordinate system (see DW97; LH07a). To obtain Qei (C; 0i 0)i Qea (C; 0i 0) ^^id (5e3(C;0i0) from Qr(6,/3, $), 
we need to use the relations between ^, -0, and 6, /?, <I> (see WD03; HL08a). 

FLIPPING PROBABILITY 

Once the mean flipping time ttf is known, one can calculate how much time the grain spends in each flipping state 
during a time interval At. 

WD03 proposed an algorithm to calculate such a flipping probability. Let Pq be the probability that the grain does 
not flip (i.e., G{t) in Appendix B), and fsame is the probability the grain stays in the same initial flipping state after 
At, they derived 

/._4(i + =!(z^^Mf!=M^)^ (Di) 

The probability that the grain in the original flipping state after At is the total probability that the grain does not 
flip plus the probability that the grain flips back original state if flip occurs 

fsa^ne = Po + (1 - Po)fl + (1 - Po)f2:; (D2) 

where is the fraction of time the grain spends in original state before the i flipping occurs. Hence, 

fsame = Pq + (1 " Po)fs (D3) 

The fraction of time At the grain spends in the initial flipping state is then 

r fsame Pq (r)A\ 

i — -TO 



where 



Po = exp(-^). (D5) 
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